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n array feed combining system for the recovery of signal-to-noise ratio (SNR) 
loss due to antenna reflector deformation has been implemented and is currently 
being evaluated on the Jet Propulsion Laboratory’s 34-meter DSS-13 antenna. In 
this system, the de focused signal Held captured by a focal plane array feed is recov- 
ered using real-time signal-processing and signal-combining techniques. The current 
signal-processing and signal-combining algorithms are optimum under the assump- 
tion that the white Gaussian noise processes in the received signals from different 
array elements are mutually uncorrelated. Experimental data at DSS 13 indicate 
that these noise processes are indeed mutually correlated. The main result of this 
article is an analytical derivation of the actual SNR performance of the current 
su optimal signal-combining algorithm in this correlated-noise environment The 
analysis here shows that the combined signal SNR can either be improved or de- 
gra ed depending on the relation between the array signal and noise correlation 
coefficient phases. Further performance improvement will require the development 
of signal-combining methods that take into account the correlated noises. 


I. Introduction 

ffr ,?r rat r ° f <ieCP_SPaCe < ,? mmUnicatl0n networks at higher carrier frequencies has the advantage of 

-uLThTt S T S 95 *1 in r ased bandwidths for enhancing telemetry capabilities. However, 

the use of higher frequencies also has certain disadvantages. These include more stringent antenna 

P mg requirements and larger receiving antenna signal-to-noise ratio (SNR) losses due to mechanical 
deformations of large reflector surfaces. These SNR losses become more significant at higher frequencies 

t h e r a C 9 IT ' ° , ''' V ' ive e r s become smaller than the mechanical imperfections of the reflector. This is 
he case in the Jet Propulsion Laboratory’s Deep Space Network plan to employ Ka-band (32-GHz) 
communications using 34- and 70-meter receiving antennas. ’ 

An array feed combining system for the recovery of the SNR loss due to antenna reflector deformation 

tlcTsedlZThelTs Aulr 6 ^ E ^ P ' ane feed arra >' is used to collect the 

etocused signal fields. All the signal power captured by the feed array is then recovered using real-time 

a'^ob^air^^ronT t hg S ^ naPc ° mbl a techniques. In phase and quadrature, baseband signal samples 
are obtained from the downconverted received signal of each of the array feed elements and then are 

recombined after application of combiner weights. The optimum combiner weights that maximize the 
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i . , ai (mnl SNR were derived in [ll under the assumption that the white Gaussian noise processes 

tigZTr « iray elements me mutually uncorrela^Theee optimum weights 
depend on unknown signal and noise parameters that need to be estimated. The work in (1) prop 
to estimate the optimum weights from the observed residual carrier received signal samples using a 
maximum likelihood (ML) estimation of these unknown parameters. The actual combined signal S 
“ulorrelated-noise environment was also derived in [1] when the estimated weights were used in 

place of the optimum weight coefficients. 

The array feed combining system is currently being evaluated at the JPL DSS-13 34-meter antenna. 
Although the work in [1] assumed mutually uncorrelated-noise processes, experimental data [2] indicat 

that the white noise processes in the received signals from different feed elements are indeed correl ^ e . 

. a* . . r of 0 01 under clear-sky conditions. Since the noise in each of 

»^o~wh. t : noise 'plus noi* due to background tadlatlom 

this small correlation is conjectured to be caused by near-field atmospheric backgroun noise. oug 
the observed correlation in [2] is quite small in the current array feed combining system, future plan 

provides an exact analysis of the combined signal SNR performance in this correlated-noise environ . 
The performance analysis here considers only the signal combining algorithm proposed in [1] which 

t^r° ef^t 1 of^ tffis^ correta^on ^is^n 1 the T£ analysis 

here shows that the combined-signal SNR can be either improved or degraded depending on the rela ion 

noise processes. Our work on this problem is still in progress. 


II. Array Feed Signals and Combining Algorithm 


’ Consider a ^-element array and the NASA Deep Space Network ataudardr^duaWmr modulation 

l£! la ,.] 

Er« 

the array elements. Let 


yk{i B ) = Videos 6 + js(i B ) sin 6} + n k (it}), ( x ) 

denote the stream B signal samples from the *th array element. The complex signal parameters 


Vfc = |Rjcl e ^ tffc 


( 2 ) 
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represent the unknown signal amplitude and phase parameters induced by the antenna reflector defor- 
mation. Moreover, <5 is the modulation index, s(i B ) = ±1 is the transmitted data, and {njt(i B )} is the 
zero-mean white Gaussian noise corruption in the stream B signal samples from the fcth array element. 
The primitive baseband signal samples are also more narrowly low-pass filtered by averaging succes- 
sive blocks of M a samples to yield a residual carrier signal stream A for each array element. Clearly 
M a > Mb, and 17 = M A /M B is the ratio of the bandwidth of stream B to stream A. Let 


“*(m) = V k cos 6 + m k (i A ), i A = 1,2, ••• (3) 

denote the stream A signal samples from the kth array element. Here {m fc (i^)} is the zero-mean white 
Gaussian noise corruption in the stream A signal samples from the fcth array element. 

Let A t and £■ denote the transpose and complex conjugate transpose of the matrix A, respectively. 
The white noise sequences corresponding to different array elements are assumed to be correlated. To 
specify these correlations, consider 


S(ifi) = (ni(i'B)r ■ ■ ,n K (i B )) T 
m{i a) = (m x {i A ),--- ,m K {i A )) T 


Then {n(i B )} and are each sequences of independent identically distributed (i.i.d.) zero-mean 

complex Gaussian random vectors of dimension K . The respective covariance matrices 


E B = { r Bkj} = E [n(fs)n(f B ) t ] 

Ea = {r A kj} = E \^{i A )m(i A )^ 

of n(i B ) and m(i A ) then specify the mutual correlations between the white noises in the signal streams 
from different array elements. For example, r Bkj is the correlation between the noise variables n k (i B ) 
and rij{i B ) in the stream B signals from the /cth and jth array elements, respectively. Moreover, define 


PBkj 


r _Bki 

\fBBkkX~Bj] 


= \PBkj\e jlpak:i 


(4) 


to be the correlation coefficient between the noise samples n k (i B ) and n 3 (i B ). We shall assume as in 
[1] that the complex Gaussian noise samples and rrik{iA) each has statistically independent real 

and imaginary parts of equal variance. This assumption is not required for the following analysis, but 
is made to maintain consistency with the results reported in [1], So, 2 a% k = r Bkk and 2 a\ k = r Akk are 
the respective variances of n k (i B ) and m k (i A ), where a\ k and a\ k are the respective variances of the 
real or imaginary parts. Because of the different averaging rates in streams A and B on the primitive 
baseband signals, it follows that R g = r/R A . Finally, these different averaging rates also imply that m{i A ) 
is independent of n(i B ) provided that i A < i B and the samples averaged to yield m(i A ) occurred prior 
to the samples averaged to yield n(i B ). 

The complex combining weight coefficients w k , 1 < jfc < K, given by 


w k 


V k * 

^\k 



(5) 
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were shown in [1] to maximize the SNR of the combiner output in the uncorrelated-noise case, resulting 
in a maximum possible SNR equal to 


= £ 


wr_ 

2 V<T 2 Ak 


V l^l 2 

h 2a b* 


( 6 ) 


That is, the optimum attainable SNR in the uncorrelated-noise case is equal to the sum of the SNRs 
of each of the feed array element outputs. The signal parameters V k and the noise variances a Ak are 
unknown parameters that need to be estimated to obtain an estimate of the optimum weight coefficients. 

Assume that these unknown parameters are not random. The estimates for V* and a Ak developed in 
[1] are univariate sampling estimates based on the stream A residual carrier signal samples {u*(m)}- In 
the uncorrelated-noise case, the stream A signal samples from different array elements are statistically 
independent. Hence, estimates of the weight coefficients w k based on these estimates of V k and a Ak are 
also mutually independent. However, in the correlated-noise environment, these signal streams are no 
longer mutually independent and, hence, the resulting estimates for w k are also no longer independent. 
In order to put this dependence in the proper perspective for the SNR performance analysis below, we 
will describe the estimation techniques developed in [1] in terms of multivariate sampling estimates based 
on the vector of stream A signal samples {u(i>i)} where 


u(i A ) = (ui(u), ■ ■ ■ ,uk{ia)) 


T 


Instead of estimating V k directly, consider estimating X k = V k cos 6. Define 

X = (X l ,---,X K ) T 


Then it follows from Eq. (3) that {u(m)} is an i.i.d. sequence of complex Gaussian random vectors with 
mean X and covariance matrix R A It follows from multivariate statistical analysis [4,5] that, based on 

observations {u(iA — !)>••• > U.{^a ~ -^)}i 


X(ia)= = ~[ Y a(J) 

1—lA ~ L 

is the ML sample mean estimate of X and 

R a {ia) = 2X^2 [- (i) _ [ m(Z) " 

l—x A -L 


(7) 


(8) 


is equal to (L - 1 )/(L - 2) times the corresponding sample covariance estimate of R A . The approach 
in [11 uses X k (i A ) as the estimate of X k and consequently V k (i A ) = X k {i A )/cos 6 as the estimate of V k . 
Moreover, the fcth diagonal element 2 & 2 Ak {i A ) of Ra^a) is used in [1] as the estimate of 2 o Ak , which is 
the kth diagonal element of R A . Finally, the estimate given by 


- r x YkM _ EM 
Wk ' lA ' 2ri& 2 Ak (i A ) 2 t] cos 6 a\ k (i a ) 


(9) 
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was shown in [1] to be an unbiased estimate of the optimum combining weight coefficient w k given by 
bq. (5) in the uncorrelated-noise case. These weight coefficient estimates are used in a sliding window 
structure to produce the following combiner output sequence: 


K 

2 (*b) = '52 Wk (m) Vk(i B ) 

k= 1 


( 10 ) 


where i A is the largest integer less than i B , so that the residual carrier signal samples 
y(i') ~ used for estimating w k (i A ) occur before the full-spectrum signal sample 


III. SNR Performance Analysis 

The objective is to determine the actual SNR of the combiner output in the correlated-noise environ- 
ment. From Eqs. (1) and (10), the combiner output can be written as 


z (iB) = s c (i B ) + n c (i B ) (H) 


where 


K 

scOb) = 52™*$ a) Vk^ 3 ^ 6 

k = i 


and 


K 

n c(*s ) = 52 ™ k (m) n fe(*B) (13) 

k = 1 


are the signal and noise components, respectively. Since the residual carrier signal samples used for the 
estimates w k (i A ) occur prior to the full spectrum signal samples y k (i B ), and since {m k {i A )} and {n,(i H U 
are l.i.d. sequences, it follows that w k (i A ) and rXj{i B ) are uncorrelated random variables for every k and 
3 ■ Each has zero mean - lt then follows from Eqs. (13) and (12) that n c (i B ) also has zero mean and 

is, moreover, uncorrelated with s c (i B ). Let Var[Z] = E [\Z - E[Z}\ 2 ] denote the variance of a complex 

random variable Z. Thus it follows from Eq. (11) that the actual SNR of the combiner signal output 
z{i B ) given by Eq. (10) can be written as 


lML = 1 E ^)1| 2 = |E[s c (f fl )]| 2 

Var[z(i B )] Var[s c (i s )] 4- Var[n c (i B )] 


(14) 


It is well known [4,5], that X_(ia) and R A (i A ) are statistically independent and that 
has a chi-square distribution with 2 (L - 1) degrees of freedom. As a result of these 
from Eq. (9) in a derivation similar to that in [1] that, for 1 < k < K, 


^L-2)a\ k {i A )la\ k 
properties, it follows 


E [^(m)] = w k 


(15) 
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where fu*} are the optimal combining weights given by Eq. (5). That is the <*timat«l weight coefficient 
are unbiased as in the uncorrelated-noise case [1]. It then follows from Eqs. (12), (15), (5), and (6) that 
for both the correlated- and uncorrelated-noise cases, 


K 


|E[s c (ifj)]| = \e js{iH)6 ^E[wk{i A )}Vk\ = 7 

k= 1 


(16) 


Consider next the variances of s c (r B ) and n c (i B ) in Eq. (14). Using Eqs. (12) and (15), we have 


= 22 Y E [OMCO - w k)(wjCi A ) - ^)*E fc U/] 


Var[s c (i B )j 


(17) 


k = 1 j=l 


where urn is given by Eq. (5). Consider first the case when the Gaussian noise processes in the signals from 
different array elements are mutually uncorrelated. Since »*(«> and *,(>„) are pairwise mdependent 
for k ^ j in this case, the variance of s c (is) can be written as 


K 

Var[/[s c (2s)] = (^)l ^ 


(18) 


Let 


■ K k 

ft = 2 He { Y, Y VkV J { E ^ ^ Ml " WkW 3 / 

k— 1 J -k ~\- 1 


(19) 


Combining Eqs. (17), (18), and (19) then yields 


Var[s c (i B )] = Var[/[s c (i B )] + ft 


( 20 ) 


Recall that n c (i B ) has zero mean and w k (i A ) is statistically independent of n,(* B ) for all k and j. Them 
similar to the derivation leading to Eq. (20), we can write 

K K 

Var[n c (i B )] = ^5;E[u) fc (i A )u)*(ft)] E [n fc ( lB )n*(i B )] = Var B [n c (z B )] + ft &1\ 

= 1 j = 1 


where 


= 2 Tie 


K 


Y Y E K (toi E n 'j 

fc=i j=fc+i 


( 22 ) 


and where 
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Var u[n c (i B )] = [| w k (i A ) f 

k = 1 


Var [|n*(i B )| 2 ] 


is the variance of n c (iB) in the uncorrelated- noise case. It then follows from Eqs. (14) and (16) that the 
actual SNR of the combiner output in the uncorrelated-noise case is given by 


Jml Var B [s c (i B )] + Var B [n c (z B )] 
So it follows from Eqs. (14), (16), (20), (21), and (23) that 

( j ^) 


(23) 


(24) 


where 


Pi 4- 02 

Var C /[s c (? J g)] + Varc/[n c (z B )] 


(25) 


The factor 1/(1 + d) in Eq. (24) represents the improvement in SNR caused by the correlation between 
the noises in the signals received from different array elements. Note in particular that 0\ and 02 can be 
either positive or negative in value. Hence, an SNR improvement is obtained when d is negative and a 
degradation is obtained otherwise. 

Expressions for Vary [s^b)] and Var^r [n c (is)] are given in [1]. Thus, we need only determine 0\ and 
02 to obtain d and thereby obtain an expression for 7 ml from Eq. (24). In order to do this, we need 
only obtain an expression for E[^(u)^(2^)] when k ^ j. Using the property that X(i A ) is statistically 
independent of R A (i A ), it then follows from Eqs. (9), (7), and (8) that, for k / j , 


E [w k (i A ) w] (i A )] = - ^ E [X’ k (i A ) Xj (u)] E 


lv 2 Ak (iA)<T 2 Aj (i A ) _ 


Since X_(i A ) has mean X_ and covariance matrix jjR A — J Z—B [5], it follows that 


E 


{*a) Xj (z^)J - — r* Bk j -f X k X j 


(26) 


(27) 


Recall that 2(L — 2)a Ak (i A ) is the kth diagonal element of the matrix ( L — 2)R A (i A ). Let 


A = 



A 12 

^22 


be a 2 x 2 matrix where An and A 22 are the kth and jth diagonal elements, respectively, and A i2 is the 
element in the fcth row and jth column of (L - 2)R A (i A ). So we have 
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(28) 


VaS a) 


= 4 (L - 2) 2 E 


1 

^11^22 


Complex multivariate statistical sampling theory [4,5] has shown that A has the same distribution as 
that of J2i=i ZlxZ-x where {ZJ is a sequence of i.i.d. zero-mean complex Gaussian random vectors with 
covariance matrix E given by 


T Akk T Akj 

1 

TBkk 

r Bkj 

_ r Akj rA Jj . 

V 

jBkj 

r Bjj _ 


(29) 


This type of distribution is called a complex Wishart distribution [4,5], with parameters E and (L - 1). 
Denote the determinant and trace of a matrix A by \A\ and tr(^4), respectively. Then if L > 4, the joint 
Wishart probability density of (-An, A 22 , -A 12 ) is given by [4] 


p(Au, A 22 , A 12 ) 


g 

(^11-^22 — |-4i2 1 2 j 
■kT{L — 1)T(L — 2)|E| t_1 


exp [— tr(E 1 ^4)] 


(30) 


for An, A 22 > 0 and \An\ 2 < A n A 2 2 , where r(x) is the gamma function. The derivation in Appendix A 
obtains the expression given by Eq. (A-5) for E[l/An A 22 ] starting from Eq. (30). Define for L > 4 and 
0 < x < 1, 


vr _, ^ fk + L- 3\ x k 
h{x) = (L - 2)(1 - x) E ( k ) k + L-2 

k = 0 v ' 


(31) 


Assume that the correlation coefficients between noise components of the fcth and jth array element 
outputs pBkj given by Eq. (4) are always less than one in magnitude. Then, by using Eqs. (28), (49), 
(31), and (27), Eq. (26) can be written as 


E [wk (m) w] (u)] = h (iPsfcjl 2 ) 


PBkj 


+ 


VW 


T]L COS 2 i5 2(J Bk& Bj ^ a Bk a Bj 


(32) 


When \p B kj\ < 1 and L > 4, we obtain, by using Eqs. (2), (4), (5), and (32) in Eqs. (19) and (22), 


0 1 + 02 


= 2 E E lfL(\pBkj\ 2 ) 

fc=lj=fc+l V 


\PBkl_ , \Vkf}Vjl 
■qL cos 2 6 4(r%k*Bj 


+ 



rjL cos 2 6 


2 CTBk&Bj 


\pBkj I COS (ti kj ~ PBkj) 


|Vtl 2 \Vj\ 2 \ 

te 2 Bk°Bj J 


(33) 


where pBkj is the phase of the correlation coefficient pBkj between n*(i B ) and rik(i B ) and where ■dkj = 
6 k - 6] is the phase difference between the signal components of the fcth and jth array elements. Finally, 
by using Eqs. (44) and (48) of [1] for Vary [s c (i B )] and Vary (n c (r B )], respectively, Eq. (25) can be written 
as 
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Pi + P2 


2 


(34) 


d = 

(r^i) h + (T + K)/riLcos 2 6\ + (jj^) 



where j3\ + P 2 is given by Eq. (33) and 7 is given by Eq. ( 6 ). In order to arrive at an explicit expression 
for 7 a fLj we note that Eqs. (44) and (48) of [ 1 ] in Eq. (23) give 


7 ML 



[7 + (y + K)/riL 



(35) 


So the actual SNR of the combiner output in the correlated-noise case can be determined from Eqs. (24), 
(34), and (35) when L > 4 and \pBkj\ < 1- The two measures of particular interest in understanding the 
SNR performance are 1/(1 + d) and 7ml / 7 * The measure 1/(1 4 * d) represents the gain in SNR caused 
by the correlation between the array element noises and will be referred to as the correlation gain. In 
the uncorrelated- noise case, 7 / 7 ml represents the loss in SNR due to the combining algorithm since 7 
is the maximum possible achievable SNR. We shall adopt the same measure here and define 7 ml/ 7 as 
the combining gain for ease of comparison with the uncorrelated-noise case. The combining gain also 
represents the gain in SNR over the sum of SNRs of the individual array element outputs. 

Let us examine the characteristics of the SNR performance. In the uncorrelated-noise case, the actual 
SNR performance 7 ^ /L converges to the maximum possible SNR achievable 7 as the number of samples 
L approaches infinity. It is interesting to also examine the combining gain in the correlated-noise case 
as the number of samples approaches infinity. It is shown in Appendix B that //,(#) — * 1 as L — > 00 for 
0 < x < 1 . Assume that the pairwise noise correlation coefficients pskj are all less than one in magnitude. 
Then, taking the limit as L — ► 00 in Eqs. (34) and (33) yields 


lim d = — 

L—*oc 7 


K K 


E E 

k = 1 1 


mvj\ 

2<^ Bk& Bj 


\PBkj\ cos {dkj ~ P Bkj ) 


(36) 


So the limiting value of d can also be of either sign, positive or negative. In fact, the limiting value is 
always negative if dkj — pBkj ~ ?r for all k j, and always positive if - (fskj = 0 for all k ^ j. 
It then follows from Eq. (24) that as L — > 00 , the limiting value of the actual SNR performance 7 ml 
in the correlated-noise case can be either greater or smaller than the maximum possible SNR 7 in the 
uncorrelated-noise case, depending on the relation between the signal and noise correlation phases. This 
is not really that surprising, since the maximum possible SNR performance in the correlated-noise case 
is generally not equal to 7 . 

Bounds on the actual SNR performance 7 a / l that depend on a fewer number of parameters than 
the exact expression are also useful. We shall derive upper and lower bounds that depend only on the 
maximum magnitude of the noise correlation coefficients and on 7 , the sum of the SNRs of the individual 
array element outputs. We first note the following inequalities derived in [ 1 ] for this purpose: 



Similar to the left-hand inequality of Eq. ( 37 ), we have 


(37) 
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Applying the left-hand inequality of Eq. (37), the inequality of Eq. (38) gets the following upper bounds. 


and 


Let 


2 


K K 


£ £ 

k— 1 j=k + 1 


\Vk\ 2 \Vj \ 2 

^ a Bk a Bj 



K K 

^£ £ 


k=lj=k+l 


JSII!5Ls(jr-ih 

Bk& Bj 


(39) 


(40) 


Pmax inaX 

be the maximum magnitude of the correlation coefficients between array element noise components. Note 
from Eq. (33) that the worst-case phase resulting in the largest possible d occurs when i) k] - = 0 

for all k ± j. Hence, application of the left-hand inequality in Eq. (37), the inequalities of Eqs. (39) and 
(40), and the bounds of Eq. (B-7) on f L {x) given in Appendix B yields the following upper bound on the 
worst-case d : 


(L - 2 )(K - 1 )jhnax [7 + {Kpmax + 7 ) cos 2 A| + 7 2 (1 - 1 /K) 
( L - 2 ) [7 + (7 + K)/rjLcos 2 6 ] + 7 2 /K 


(41) 


Similarly, since the best-case phase resulting in the most negative possible d occurs when d kj - <p B kj = * 
for all A: ^ j, the following lower bound on the best-case d can be obtained: 


(L - 2 ){K - l]Prnaxl [l + l/^LcOS 2 ^] 
“ (L - 2) [7 + (7 + K)/rjLcos 2 6] + j 2 /K 


(42) 


Finally, using the inequalities of Eq. (37) in Eq. (35) yields the following bounds on the actual SNR 
performance 7 ^ in the uncorrelated-noise case: 


7 ml 


< 


(£- 3 ) 7 2 

(L - 2) [7 + (7 + K)/r]L cos 2 « 5 ] + 7 2 /K 


(43) 


and 


jj 

1ml 


> 


(L - 3 ) 7 2 

( L - 2) [7 + (7 + K)jr\L cos 2 6 ] + 7 2 


(44) 


An upper bound on the actual SNR performance 7 a/ i. is obtained by using the lower bound of Eq. (42) 
on d and the upper bound of Eq. (43) on 7 ^ in Eq. (24). Similarly, a lower bound on 7 ml is obtained 
by using, instead, the upper bound of Eq. (41) on d and the lower bound of Eq. (44) on 7 ^^. 



IV. Numerical Example 

We consider here the numerical example in [1] of using a K = 7 element array feed in the JPL 
Deep Space Network. In this example, a modulation index 6 = 80 deg and a primitive sample period 
Tq = 2.5 x 10 -8 s are assumed. The full-spectrum modulation signal is assumed to be of bandwidth 2 x 10 6 
Hz, which yields Mq — 20. Moreover, the ratio of the full-spectrum bandwidth to the residual carrier 
bandwidth T} = M a /Mq = 200. Nominal Pt/Nq of 55 and 65 dB-Hz are considered with corresponding 
7 = {Pt/No)MbTq. Upper and lower bounds on the combining gain 7 ml/t are shown in Fig. 1 as a 
function of the number of samples L averaged to obtain the weight estimates. Here Pt/Nq = 55 dB-Hz, 
and maximum correlation coefficient magnitudes p max of 0.01 and 0.02 are considered. Convergence 
of these bounds to within 0.01 dB of their limiting values occurs at about L = 3000 samples. This 
corresponds to an averaging time of M A T 0 L = 0.3 s and supports real-time operations for antenna 
deformation compensation. The limiting upper bounds on the combining gain are about 0.26 and 0.56 dB 
for Pmax equal to 0.01 and 0.02, respectively. The corresponding lower bounds on the combining gain 
are —0.26 and —0.50 dB, respectively. The actual limiting value for the combining gain, which is given 
by Eq. (36), will fall between these bounds. Similar results are shown in Fig. 2 for Pt/Nq = 65 dB-Hz, 
where convergence of the bounds occurs at smaller values of L to virtually the same limiting values as 
the Pt/Nq = 55 dB-Hz case. 



Fig. 1. Combining gain versus L for P^Nq = 55 dB-Hz. 

Figures 3 and 4 plot upper and lower bounds on the correlation gain 1/(1 + d) for p max equal to 
0.01 and 0.02. Figure 3 considers Pt/Nq = 55 dB-Hz and Fig. 4 considers Pt/Nq = 65 dB-Hz. The 
limiting values of these bounds are identical to the limiting values of the corresponding bounds on the 
combining gain. The differences between the behavior of the lower bounds at Pt/Nq — 55 dB-Hz and 
those at Pt/Nq = 65 dB-Hz are due to the looseness of these lower bounds at small values of L. For a 
large number L of samples, the upper and lower bounds on the combining gain diverge as the maximum 
correlation coefficient magnitude increases. This can be seen from Fig. 5, which shows the upper and 
lower bounds on combining gain for Pt/Nq = 55 dB-Hz at L = 5000 samples as p max increases from 
0.01 to 0.1. The upper bound increases from 0.26 to 3.96 dB and the lower bound decreases from —0.26 
to -2.05 dB in this range of Pmax ■ The observed correlation coefficients of 0.01 magnitude in [2] were 
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Fig. 4. Correlation gain versus L for PjJNq = 65 dB-Hz. 



Fig. 5. Combining gain versus p majf for PjtN Q = 55 dB-Hz and L = 5000. 
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obtained in clear-sky conditions with a receiver noise temperature of 90 K and a system noise temperature 
of 120 K. Improvement of the receiver noise temperature to 25 K will increase the correlation coefficient 
magnitude to about 0.02. As noted above, a maximum possible improvement of 0.56 dB and a maximum 
possible degradation of -0.50 dB results. Preliminary measurements at DSS 13 indicate that even larger 
amounts of correlation occur under adverse weather conditions. This will result in even larger potentia 
improvement or degradation of SNR performance relative to the uncorrelated-noise case. 


V. Conclusion 

An array feed combiner system for the recovery of SNR loss due to antenna reflector deformation has 
been implemented and is currently being evaluated on the Jet Propulsion Laboratory 34-metei • DSS-13 
antenna. The current signal-combining algorithms are optimum under the assumption that the white 
Gaussian noise processes in the received signals from different array elements are uncorrelated. Exper- 
imental data at DSS 13 indicate that these noise processes are indeed mutually correlated. The m 
result of this article is an analytical derivation of the actual SNR performance of the current subopti- 
mal signal-combining algorithm in this correlated-noise environment. The analysis here shows that the 
combined-signal SNR can be either improved or degraded depending on the relation between the array 
signal and noise correlation coefficient phases. Further performance improvement will require the devel- 
opment of effective combining systems that take into account the correlations between the array feed 

element noise processes. 
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Appendix A 


Derivation of E 


i 

^11^22 


We first obtain the joint probability density function p{An,A 22 ) of (An, A 22 ) by integrating Eq. (30) 
over the complex region 5 = {A l2 : |Ai 2 | < \A4 11 ^ 22 } of values taken on by A l2 . Let G = {G tJ } = E" 1 
and convert the variables G 12 and A X2 into polar coordinates: G l2 = |G'i 2 |e-^ ^ and A 12 = re j<l> . Then it 
follows from Eq. (30) that 


P{A 11 , ^ 22 ) = 


0 — (Gi \A 11 +^ 22 ^ 22 ) 
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r(L-i)r(L-2)|E| i-1 


\A4i 1 A22 

J r(A ll A 22 -r 2 ) L ~ 3 h(2r\G 12 })dr 
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(A-l) 


where Iq{x) is the zero-order modified Bessel function of the first kind, which has series representation 


Io(x) = 

k = 0 


-2fc 


2 2k (k\y 


(A-2) 


By making a change of the variable of integration using the series of Eq. (A-2) and the integral relation 
(3.251) of [6], the integral in Eq. (A-l) can be written as 


y/Aii a 22 

j r(A u A 22 - r 2 ) 1 ” 3 1 0 (2r\G u \)dr = (A n A 22 ) L ~ 2 f^ 
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(A-3) 


Substituting Eq. (A-3) into Eq. (A-l) and using the fact that T(n) = (n - 1)! for integer n, we obtain 

p(A n ,A*) = Mu f MnjalgajY (A . 4) 

(t-2)!|SI t!(i + 1-2)! 1 1 
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Using Eq. (A-4), integrating term by term in the series, and obtaining G — E 1 and |E| directly from 
Eq. (29) in terms of p B kj, r B kk, and r B jj . we have 


E 


1 

AhA 2 2 


1 00 

l y 

(L-2)\n L ~ l {G u G22) L ~ 2 fo 




k + L-2 


r? 2 (l - \fBkjf) f' ^ _ 3 \ 

(L - 2)r B kkVBjj j“ V k ) k + L-2 


(A-5) 


The series in Eq. (A-5) can be shown to converge by using the ratio convergence test whenever \pBkj\ < 1 


Appendix B 
Bounds on f L (x) 


Let L > 4 and 0 < x < 1. We will obtain upper and lower bounds on /l(z) that are asymptotically 
tight in the limit as L — ► oo. First note that 

L- 2 _( L-3 \ (L- 2\ (k + L-3\ ( B -l) 

k + L-2 \k + L-3) \L-Z) \k + L-2) 


and that for A; > 0, 


£~ 3 < fc + L ~ 3 < i (B-2) 

L - 2 - k+L-2 ~ 


Using these bounds of Eq. (B-2) in Eq. (B-l), we have 

£ ~ 3 < L ~ 2 < ( L ~ 2 \ L ~ 3 _ (B-3) 

k + L- 3 “ k + L-2 - \L — 3J k + L — 3 


Next, by using the bounds of Eq. (B-3) in Eq. (31), we get the following. 
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It can be shown in [7] that for 0 < x < 1, 


oo 


E 


k + L - 4 
k 



L - 3 


Using Eq. (B-6) in Eqs. (B-4) and (B-5), we then obtain, for 0 < x < 1, 


1 < h(x) < 


L - 2 
L - 3 


The upper and lower bounds given in Eq. (B-7) are both asymptotically tight in the limit as L 
we can conclude that for 0 < x < 1, /l(^) — * 1 as L — » oo, where the convergence is uniform in 


(B-6) 


(B-7) 
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We consider the problem of finding a trellis for a linear block code that minimizes 
one or more measures of trellis complexity for a fixed permutation of the code. We 
examine constraints on trellises, including relationships between the minima] trellis 
of a code and that of the dual code. We identify the primitive structures that can 
appear in a minima] trellis and relate this to those for the minima] trellis of the 
dual code. 


I. Introduction 

Every linear block code can be represented by a minimal trellis, originally introduced by Bahl et al. [1], 
which is a labeled graph that can be used as a template for encoding or decoding. As shown by McEliece, 1 
the minimal trellis simultaneously minimizes the maximum number of states, the total numbers of vertices 
and edges in the trellis, and the total numbers of additions and path comparisons required for decoding 
with the Viterbi algorithm. 

In this article, we examine properties of the minimal trellis representation of a code and its dual for 
a fixed permutation. A companion article [2] uses these results to examine the problem of finding a 
permutation that minimizes one or more trellis complexity measures. 


Section II reviews the subject of minimal trellises for a fixed permutation of a code. We examine the 
building blocks of such trellises and identify several different measures of trellis size or complexity. In 
Section III, we illustrate the connection between the minimal trellis of a code and that of the dual code. 
The section includes results that describe the structure and complexity of trellises for self-dual and other 
special codes. 


1 R. J. McEliece, “On The BCJR Trellis for Linear Block Codes,” submitted to IEEE Trans. Inform. Theory. 
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